Universal Properties of the Ultra-Cold Fermi Gas 
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We present some general considerations on the properties of a two-component ultra-cold Fermi 
gas along the BEC-BCS crossover. It is shown that the interaction energy and the ground state 
energy can be written in terms of a single dimensionless function /i(^,t), where ^ = —{kpas)'^ 
and T = T/Tf- The function h{^,T) incorporates all the many-body physics and naturally occurs 
in other physical quantities as well. In particular, we show that the RF-spectroscopy shift 5uj{£,,t) 
and the molecular fraction /c(^,t) in the closed channel can be expressed in terms of /i(^,r) and 
thus have identical temperature dependence. The conclusions should have testable consequences in 
future experiments. 

PACS numbers: 



I. INTRODUCTION 

Over the past few decades, there have been consider- 
able efforts and progress in understanding the physics 
of ultra-cold Fermi gases[l-3, 6]. In general, the the- 
oretical investigations fall into two categories depend- 
ing on how one incorporates the physics of Feshbach 
resonances. In the so-called single-channel model, one 
neglects the closed channel component, while incorpo- 
rating its effects through the open channel scattering 
length a^, given by. 



as{B) = flfcg 1 



AB 

B-Bn 



(1) 



where a^g is the background scattering length in the 
absence of inter-channel coupling and B is the exter- 
nal magnetic field. AB is the width of the resonance 
and Bo is the position of the resonance. See Table. V. 
in Ref.[5] for a list of values of these parameters for 
the alkali elements currently under investigation. The 
approximation is valid in the case of a so-called broad 
resonance where ep *C Sc- Here is the Fermi energy 
of the system and Sc = 2^2/^^l [7] , characterizing the 

typical energy scale associated with the two-body Fesh- 
bach resonance [12]. A/i is the magnetic moment differ- 
ence between open and closed channels. The physics of 
the single channel model is essentially the same as in the 
crossover model studied decades ago in the literature [8- 
10]. Most experimental systems {e.g. ^Li at magnetic 
field B — 8340) fall into this category. In general, any 
system with sufficiently low density will be described by 
a single channel model. 

The problem associated with single-channel model is 
easily stated: Given a spin-1/2 Fermi gas with both 
spin components equally populated, with interactions 
only between opposite spin states, characterized by the 
s-wave scattering length as{B), what are the ground 
state and thermodynamic properties? At zero temper- 
ature, neglecting finite range corrections, the only rele- 
vant parameter is ^ = —{kFas)~^ . The basic question 
concerning the static properties of the system is how to 
find the ground state energy E{^). At finite tempera- 
ture, the most important questions are the calculation 
of the thermodynamic potential and in particular the 



location of the phase transition boundary Tc{S,) as a 
function of see Ref. [11] and references therein. Up 
to now, these problems have not been amenable to an- 
alytic solutions; our most reliable knowledge of those 
quantities comes from Monte Carlo simulations [13-19]. 

For most experiments, it is sufficient to use the single- 
channel model as the theoretical framework to interpret 
them. However, there are certain cases where we are in- 
terest in the physics of the closed channel specifically for 
a broad Feshbach resonance, as in the case of the Rice 
experiment to be discussed later[39]. Furthermore, in 
the so-called narrow resonance case where ep ^ ^c, the 
molecular states in the closed channel cannot be ne- 
glected. For these reasons, we have to start with the 
more general two-channel model[4, 20, 21], which in- 
corporates the closed channel on the same footing as 
the open channel. Typically, one introduces a bosonic 
operator 'l'(r), which creates a closed channel molecule 
with center of mass position r and incorporates its ef- 
fects through the coupling term in the Hamiltonian, 

(2) 



~g{^\r)^^ir)^^{r)+R.C. 



where g is the bare coupling constant. Note that the 
coupling scheme above enforces the momentum conser- 
vation in the conversion processes and thus the mo- 
mentum distribution of closed channel molecules is inti- 
mately connected with the open channel pair states. We 
shall return to this point later. Note also that the in- 
ternal structure of the closed channel molecule is frozen 
as a result of its high internal excitation energy which 
is much larger than any other energy scales relevant for 
the many-body physics. 

In general, exact solutions or even general statements 
about the above two models are quite difficult. How- 
ever, at resonance, i.e. ^ = 0, as drops out of the 
problem and we are left with only two energy scales, 
ep = S^/2m(37r^n)^/'^ and the temperature T. Here 
m is the mass of the atom and n is the density of the 
system. Thus, it follows from dimensional analysis that 
one may write the average single particle energy of the 
system at finite temperature as = 0,t) = epfEU = 
0, r), where fsi^ = 0, r) is a dimensionless function and 
T = T/Tp = ksT/ep, where fcs is the Boltzmann con- 
stant. In particular, at r = 0, the average single particle 
energy is proportional to cp, with universal constant 
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/^{^ = 0,T = 0) = 1(1 + /3). The parameter /3 has 
been calculated in many ways in the literature and is in 
good agreement with experiments. See Table. II. in [2] 
for a summary of values of /3 obtained theoretically and 
experimentally. By the same argument one can write 
down other thermodynamic quantities of the system, 
with the conclusion that the thermodynamic properties 
of the system are universal regardless of the particular 
system under investigation [23]. 

The universal thermodynamics works only at unitar- 
ity. However, it is possible to generalize the idea to 
the parameter space where ^ 7^ 0. It is recognized 
that for d > 2, the system is controlled by an unstable 
fixed point which resides near the Feshbach resonance 
in the zero density limit with attractive interactions [24]. 
By utilizing a large-N expansion as applied to Sp(2N) 
model, one can calculate the scaling form of the canon- 
ical free energy of the system. The e-expansion has 
also been used to investigate the properties of the two- 
component Fermi gas away from rcsonance[28, 29]. 

There is however, another form of 'universality' which 
is more deeply rooted in the actual physical properties 
of the system. For example, by examining the short- 
range form of the many-body wave function, one can 
show that several physical quantities (interaction energy, 
RF-spectroscopy shift etc.) depend on temperature r 
through one universal function h{^,T). This univer- 
sal dependence comes about because of one peculiar 
property of the dilute Fermi gases: the range of the 
interaction is much smaller than the inter-particle dis- 
tance, i.e. kpTo ^ 1, where rg is the range of the 
potential. Thus, important effects associated with in- 
teractions come mostly from two-body encounters. The 
argument presented below can then be regarded as an 
expansion in terms of kpTo- Let us note that this argu- 
ment can be trivially modified in the case of the imbal- 
anced Fermi gas. 

The organization of the paper is the following. In 
Sec. II, we give a general discussion of the physical sys- 
tem in terms of the two-body density matrix and sep- 
arate the two-body and many-body contributions in 
it. In Sec. Ill, we apply the rcsidt of Sec. II to several 
physical quantities and show that they can be written 
in terms of one universal function /i(^, r) which car- 
ries all the many-body dependence. There are resid- 
ual ^-dependences as a result of the two-body physics, 
which can in principle be calculated without any refer- 
ence to the many-body system. The temperature de- 
pendence of those physical quantities is universal and 
has experimental consequences as described in Sec. III. 
In Sec. IV, the main conclusions of the paper are sum- 
marized and discussed. In the Appendix, we give an 
alternative derivation of the linear ^-dependence of the 
energy of the system away from resonance, based on the 
many-body wave functions. 



II. GENERAL SETUP 

The difficulties involved in analyzing either the single- 
channel or the two-channel model are often expressed 



as a lack of small parameter because of the resonant 
interaction condition na^ 3> 1, which prevents a rel- 
atively straight-forward perturbation calculation as in 
the classic dilute Fermi gas[25-27]. The challenge lies 
in the correct implementation of the two-body physics, 
characterized by the 'dangerous' diverging scattering 
length Us, into the many-body calculations. One way 
to circumvent the difficulty is to devise other small pa- 
rameter, as in the e-cxpansion[28] or l/A^-expansion[24] 
utilized in recent work. On the other hand, even 
though the parameter na^ 1, we still have the 
small parameter kpro, where tq is the range of the 
potential. Note that in most investigations using the 
single- or two-channel model, the zero-range limit has 
already been taken; an exception is the case of a nar- 
row resonance [22]. In the following, we will try to set up 
an approximation scheme which utilizes the smallness 
oi kprQ. Even though it does not yield immediately a 
computational tool for the values of specific constants, 
say /3, it does lead to some general conclusions indepen- 
dent of the approximation scheme employed in a specific 
investigation. 

To motivate our discussion, let us consider the many- 
body wave function for a spin- 1/2 Fermi system with 
both spin components equally populated. We denote 
the total number of atoms N . For example, we consider 
a collection of ®Li atoms in their lowest two hyperfine- 
Zeeman states(|l) and |2)). Let us write down its many- 
body wave function as ^'(riCTi, r2C2 ■ ■ ■ ^^(^n)- Now, we 
separate two atoms, say atom 1 and atom 2, with op- 
posite spin orientations, far (compared to ro) from all 
the other atoms (3, 4, • • • , iV), and ask what is the form 
of the many-body wave function, when the distance be- 
tween ri and r2 is taken to lie within the range of two- 
body interacting potential, i.e., [ri — '^I'o- Since all 
the other N — 2 atoms cease to interact with atom 1 
and atom 2, we conclude that 

lim *(ri<Ti, f2a-2 • • • fjva-jv) (3) 

\ri-r2\<ro 

OC A(l){fi - r2)Si2*'(f3Cr3, ^40-4 • • • f]v0-jv) 

Here A is the trivial anti-symmetrization operator and 
E12 is the spin wave function of atoms 1 and 2. <;i(fi — 
7^2) is determined by the two-body interaction potential 
in the range jri — r2| < tq. In order for the other N — 2 
atoms to affect the form of (j){ri — f^), it is necessary 
that a third atom is at a distance < Tq. Such process 
is highly unlikely for two reasons. Firstly, the phase 
space of such event is down by at least a factor (fc^ro)^ 
as compared with the two-body encounters. Secondly, 
in a spin-1/2 system, two of the three atoms close to- 
gether miist have the same spin, thus Pauli exclusion 
principle will prevent such process from occurring. In 
fact, to the extent that we can work entirely in terms 
of the s-wavc scattering length a^, we have automat- 
ically neglected the contributions from higher partial 
wave scattering, which are of order {kproY or higher. 
We thus conclude that, to order kpro, the short-range 
behavior of the many-body wave function(in particular 
its nodal structure) is determined by two-body physics. 
An important question pertains to the form of two-body 
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wave function (l){fi — fz), since in general, the two-body 
potential can host several bound states. Here we note 
that since the many-body energy scale is much smaller 
than the energy splitting between the different energy 
levels in the potential well, it is easy to convince oneself 
that only the bound state that is closest to the scat- 
tering continuum is relevant. This state is nothing but 
the molecular state on the BEC side of the crossover. 
Here we emphasize that even in the BCS-side, where 
the two-body bound state emerges above the scattering 
continuum, the above conclusion still holds. To con- 
clude this intuitive discussion, we must point out that 
at this stage, we do not yet know the normalization of 
the short-range part of the many-body wave functions. 
This will be determined by the many-body physics. 

To make the above argument more precise, let us con- 
sider the two-body density matrix for a generic many- 
body system. The definition of the two-body density 
matrix is given by[12] 

p{fia,f2(3;f3'y,f4S;t) (4) 

Here Va(fi)i) is the Heisenberg field operator for a 
fermion with spin a. In the following, we shall con- 
sider only an equilibrium situation and thus drop the 
time t from the above expression. By the Hermicity 
property of the density matrix, we can decompose the 
two-body density matrix in the following form [12], 

/9(fiQ:,f2/3;f37,f4(5) (5) 

i 

The eigenvalues rii and eigenfunctions 0^*^(fi,f2) sat- 
isfy the following conditions, J2i '^i = N{N — 1) and 

Eo.0ld'rild'r24l*ir2,n)<t>fjrur2) = S,j. As dis- 
cussed above, in the case of a dilute Fermi gas, the 
only relevant parameter is ^ = —{kpcis)'^ and t. Thus 
Ui and (j)i3a{r2,fi) will depend on ^ and r parametri- 
cally. Now, by the argument given above in terms of the 
many-body wave function, we see that the short-range 
form of (j>ai3{ri — f2) will be determined by two-body 
physics, while many-body physics will determine the 
eigenvalues rii and the long range part of the eigenfunc- 
tions. Our philosophy in the following will be to express 
several physical quantities in terms of the two-body den- 
sity matrix and use the above facts to extract their 
universal dependence on temperature. To be success- 
ful, we need our expressions to pick up only the short 
range part of the two-body density matrix so that all 
the temperature dependence will be carried by n^'s and 
the normalizations for the pair wave functions. Physi- 
cally, as we change the temperature and the interaction 
strength, the occupation numbers of the pair wave func- 
tions (/)^*^ change while the short range part of the pair 
wave function remains the same. Let us thus consider 
an arbitrary short-range(~ ro) function s(ri — r2) and 



consider the integral, 

j dfldf2S{fl - r2) {^V'l(fl)V'2(^2)V'2(f2)l/'l(n)^ 
= Y'Tlii^^T) I S(n -f2)|^i2(fl,f2)|^dfidf2. (6) 

It is clear that in the above equation, we need only 
retain the s-wave part of the pair wave function, 
since higher partial waves have vanishing probabil- 
ity at the origin and thus hardly contribute to the 
above integral. Thus, we can write ^i2(ri,r2) = 
^-i/2gip--L_-^^W^^-)Y-^^^^^ .^vJicrc r = |fi - and Yi^ 
is the Z = m = spherical harmonics. The factor 
Q-i/2gip. 1^ ^ describes the center of mass motion of 
the pair state. Then the R.H.S. of Eqn.(6) can be writ- 
ten as, 

^n,(C,T) / drs{r)\x^l{r)\\ (7) 

By the above argument, Xi2(r) will have the form of 
the two-body radial wave function at short distance. In 
particular, in the region where as,kp^ r > ro, we 

can write X12W = C^''\^,T)xi2{r), where, 

X«(r)^l-f (8) 

has been normalized in such a way that it approaches 
1 in the region as,kp^ » r > ro. The C^'^'s are some 
constants which in principle depend on the many-body 
physics. Eqn.(6) can then be written as, 

J dfidf2S{fi - r2) ^Vl(fl)V'2(^2)V'2(f2)V'l(fl)) 

= ^nW(C,T)|C«(^,r)|2y drs{r)\x,2{r)\' 

= h{^,T)kFN I drs{r)\xi2{r)\\ (9) 

where we have defined a positive definite universal func- 
tion 

h{^,T)^Y.^^^^Mr\^^''(^'^^\">^- ^^^^ 

i 

The factor kp is inserted in order to make h{^, t) a di- 
mcnsionless function. The integral in Eqn.(9) is a con- 
stant depending on the function s(r), but it is purely a 
two-body quantity and can be calculated without mak- 
ing reference to the many-body system. In particular, 
it does not depend on the temperature T. We further 
note that the integral displays no singular dependence 
on Gs as we approaches the resonance. Thus for the 
discussion of many-body physics, it can be regarded as 
a known parameter. The intricate many-body correla- 
tions are then incorporated in one universal function 
h{^, r) and are themselves universal. As we shall show 
later, at unitarity, h{£, = 0,r) must be finite and thus 
we conclude Eq.(9) scales with kp &t unitarity. 
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Before ending the discussion of this section, let us re- 
mind ourselves of the assumptions made so far: 

a), Only s-wave scattering is important. The neglect 
of higher angular momentum(/iZ) partial waves is jus- 
tified because there are of relative order {kpro)^'' and 
thus negligible as compared with s-wave scattering. In 
fact, in the model Hamiltonian considered in the liter- 
ature, only s-wave scattering is included. 

P), The short-range form of the pair function x{^i ^ 
is determined by two-body physics and moreover, 
corresponds to only one particular form of the two-body 
wave function in the range ^ tq. The former assump- 
tion is justified because of 1), the diluteness of the sys- 
tem kpro ^ 1 and 2), the 'exchange hole': the Pauli 
principle forbids two particles with like spin to be close 
to each other. The later assumption come from ener- 
getic considerations: as long as we are interested in the 
many-body physics, which has a typical energy scale e^, 
the relevant two-body state is the one that is closest to 
the zero-energy scattering state, with all the other two- 
body states far away to be of any practical importance. 

In the following, we shall consider a uniform system 
with density n at temperature T. The interactions be- 
tween particles can be written as 

I E (/(^"^ - + 9iri - fj)Si . Sj) . (11) 

Here /(r) and g{r) are the direct and exchange interac- 
tion respectively. S is the spin operator of the valence 
electron of the atom under consideration. Experimen- 
tally, one normally works with an equal population of 
atoms(say ^Li) in the lowest two hyperfine states 1 1) and 
1 2). To the extent that one can neglect the closed chan- 
nel component, as is the case for a broad resonance, one 
may replace the full interaction by an effective short- 
range interaction in the open channel V\{r) = V{r,X), 
where A is a controlling parameter by which one can 
tune the scattering length as[12]. However, in dis- 
cussing the closed channel population, it is necessary 
to introduce explicitly the inter-channel coupling term 
W{r) which converts open channel pair states to closed 
channel molecules, see the discussion in Sec. III. D. 



III. PHYSICAL QUANTITIES 

In the following we shall discuss several physical 
quantities that can be expressed in terms of the 
universal function ft (^, r) and thus display a universal 
dependence on the temperature T. 



A. Interaction Energy 

The simplest physical quantity that can be cast in the 
form of Eq.(9) is the interaction energy of the system. 
According to the discussion above, since the interac- 
tions between particles are of short-range form, we can 



write the interaction energy per particle as, 

J dfidf2y(fi - f2) ^V'I(?l)V'2(^2)V'2(r'2)V'l(n) 

= Cv{as)kFh{tT), (12) 

where Cy(as) = / drV {r)\xi2{'>')\^ is a well-defined, 
purely two-body quantity. All the many-body depen- 
dence of interaction energy is encapsulated in the uni- 
versal function h{£,,T). Since the interaction energy 
must be well-defined at imitarity, /i(^, r) must be finite, 
and moreover free of any divergence as ^ approaches 
zero. Thus the average interaction energy per parti- 
cle scales as kp at unitarity. The interaction en- 
ergy of the system depends on microscopic details of 
the system, even at unitarity, as is clear from the factor 
Cv(as)- This result should be compared with the total 
energy of the system, to be discussed in the next sub- 
section, which is proportional to the Fermi energy ep 
at unitarity, independent of microscopic details. 



B. Total Energy 

To derive an expression for the total energy of the 
system, we first recall that, if A is the tuning parameter 
of the potential by which the scattering length can be 
varied (see Appendix), then the relation between and 
A is given by, 



Sa, 



' dr^^^lxMl' ] 6X. (13) 



On the other hand, according to the Hellmann-Feynman 
theorem, we have. 



_ / dVjr, A) 

ax ~ \ dx 



(14) 



with the average taken over the many-body state as in 

av 

OA 



(6). Since ^ is a short-ranged function, we can use 
Eqn.(13) to rewrite Eqn.(14) in terms of aj^; we find, 

(15) 



a p fc2 

J=--Nkph{^,T) 

J- fim 



da. 



m 



Here we have used the definition of /i(^,t) in Eq.(9). 
Or in terms of ^ = — (/cfOs)^^, 



dE 



m 



(16) 



Since by definition h{^,T) is a positive definite function, 
we find the somewhat trivial result that the ground 
state energy is a monotonically increasing hmction of 
^. The boundary condition on the above differential 
equation is easily obtained. Consider the case when 
^ — +00 and t — 0, we then have a free Fermi gas with 
the average single particle energy e(^ = oo) = ^ = |eir- 
Integrating Eqn.(16), we find that the single particle en- 
ergy at zero temperature along the BEC-BCS crossover 
is given by 



<0 



:£F — 2eF 



r 



(17) 
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where e(C) = e(C,r = 0) and h{£,) = h{S,,T = 0). In- 
tuitively, h{^) accounts for the reduction of the single 
particle energy due to interaction effects. At unitarity 
eo = = 0,r = 0) = (1 + /3)|eF, so we find 
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(18) 



The generalization of the above expression to finite 
temperature is straightforward but may be less useful. 
Around unitarity where ^ <C; 1, we can obtain an ex- 
pansion of the average single particle energy at finite 
temperature in term of ^. To this end, we can integrate 
Eq.(16) from ^' = to ^' = ^ and we find. 



r) - = 0, r) = 2eF t h{^' , T)d^' 
Jo 



(19) 



For ^ close to zero, the question reduces to the expan- 
sion of h{^). From the discussion in the previous sub- 
section, we know h{(, = 0,t) is finite, so we conclude 
the energy correction away from unitarity is linear in ^ 

and given by. 



n is kept constant. In the extreme BEC limit, the single 
particle energy equals the chemical potential: e = n = 
~2maJ- verifies that this is satisfied by Eq.(22). In 
the BCS limit, the first order correction to the energy 
will be of order as, namely of order coming from 
Hartree-Fock corrections. Since = — ^~^g|rT = 

—as-§^ ^ when as 0, we thus recover the usual 

relation between chemical potential /x and the average 
single particle energy of the free Fermi gas e = At 
unitarity, ^ = 0, we find again the free Fermi gas result 
e = |/U, if we assume that the energy is continuous at 
unitarity. 

Finally, using Eqn.(16), we can write the zero tem- 
perature chemical potential in terms of h{$) as, 

m(0 = (l - \m) hiOdi'^ . (23) 



Setting ^ = 0, we recover Eq.(18) since at unitarity 
/x= (l + /3)eF. 



e(C, t) - eiC = 0, r) = 2eFh{^ = 0, t)C + ■ 



(20) 



At zero temperature, the value of = 0, r = 0) can 
be calculated using the e-expansion, where one has to 
sum over all the higher order logarithms in order to re- 
cover the correct linear ^-dependence of the energy [29]. 
Instead of referring back to the conditions imposed by 
the interaction energy on the function /i(^,r), one can 
have a direct derivation of the linear ^-dependence of 
the energy away from unitarity by a straightforward 
generalization of the argument in the two-body case. 
This is presented in the Appendix. 

Before we conclude this subsection, we would like to 
derive a simple relation between the chemical poten- 
tial 12 and the average single particle energy e at zero 
temperature and thus enable us to write down the zero 
temperature chemical potential in terms of h{^). Note 
that at T = T = 0, we can write the single particle en- 
ergy as e = spfEiO- Thus using the thermodynamic 
relation P = and E = Ne, we find. 



p = n — 
on 



dn 



(21) 



f.3 

Using the relation n = we can make a change of 
variable to kp and write Eq.(21) asp = ^nkp-^^- Now, 
using the expression e = epfiO ^^'^ the fact that /(^) 
only depends on the combination ^ = — (fci^a,,)"^, we 
obtain, p = |ne -|- ^nas^p. We consider the situation 
when the density of the system is fixed and write the 
above expression as e = + ^^f|. At T = 0, we have 
the thermodynamic relation, e = where p is the 

pressure and n is the average density. We find, 
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(22) 



The above expression is very general and works along 
the whole BEC-BCS crossover provided that the density 



C. RF-spectroscopy Shift Suj 

One of the early experiments which indicated the ap- 
pearance of a new low temperature quantum state in 
ultra-cold Fermi gases was the radio-frequency spec- 
troscopic experiment carried out by the Innsbruck 
Group[30]. The experiment works with the lowest two 
hyperfine-Zeeman states(|l) and |2)) of ^Li. A radio- 
frequency field is applied to drive atoms from state |2) to 
|3). It is found that at high temperature, the frequency 
of the rf-field coincides with the bare atomic transition 
from |2) to |3), while at low temperature, there is an up- 
shift in the rf- frequency which indicates that the system 
is in a new quantum state. It is now understood that 
the full understanding of the RF-spectroscopic profile is 
quite complicated, requiring a proper treatment of the 
final-state interactions [3 1-38]. It has been shown that 
the average shift in the RF-spectroscopy is given by the 
following expression[38], 

^ GiH) + J{H) 



5U! 



hN2 

X j 5(n -r2)<^Vi(^i)V'2(^2)V'2(r2)V'i(r'i) 
G{H) + J{H) 



h 



-Cgias)2kFh{i,T) 



(24) 



where Cg{as) = J g{T)\xi-2{''')\^ dr and the functions 
G{H) and J{H) are given in Ref. [38], g{r) is the ex- 
change interaction and N2 = A^/2 is the particle number 
in hyperfine-Zeeman state |2). Again Cg{as) is indepen- 
dent of temperature T and we conclude that the average 
RF-shift has the same temperature dependence as the 
interaction energy at arbitrary ^. Note that 5u) scales 
with kp unitarity. 
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D. Closed Channel Fraction 

One of the key physical quantities in the BEC-BCS 

crossover using Fcshbach Resonance is the population 
in the closed channel. This quantity has been experi- 
mentally determined using an optical molecular spectro- 
scopic technique [39]. The experiment uses ^Li atoms in 
their lowest two hyperfine states |1) and |2), in which 
they interact primarily through the electronic triplet 
potential. On the other hand, the associated closed 
channel molecules induced by the Feshbach resonance 
interact primarily through the much deeper electronic 
singlet potential. In the experiment, a laser beam in- 
duces an electric dipole transition between the closed 
channel molecular state(X^E+, = 38) to another 
closed channel molecular state with v = 68,^^1]+. 
At low temperature, it is inferred from the loss signal 
that on the BCS side of the resonance, there is a fi- 
nite fraction of closed channel molecules which is not 
supported by the two-body physics. Thus it is sug- 
gested that the many-body quantum state must have 
non-trivial two-particle correlations like those in the 



J 



BCS state to account for the observed one-body de- 
cay in the BCS side[39]. To address the closed chan- 
nel fraction theoretically [40-42], let us first identify the 
inter-channel coupling W{r) from the bare interactions 
between the two atoms, U{r) = f{'r)+g{fi)Si ■S2- f{r) 
and g{r) are the direct and exchange interaction respec- 
tively. We restrict ourselves to the case when there are 
only two channels involved, namely, an open channel 
with atoms in the lowest two hyperfine-Zeeman states 
|1) and 1 2) and the corresponding closed channel with 
atoms in hyperfine-Zeeman states [1) and |4). Notice 
that one of the hyperfine-Zeeman states is common to 
the open and closed channels. We shall denote the in- 
teraction potential in the open and closed channel by 
Vo{r) = (12|t/|12) and K(rO = (14|?7|14) respectively. 
Now, the inter-channel coupling can be written as. 



W{n - r2) = gin - f2)(l4|5i • S2I12), 



(25) 



where |a/3) denotes a spin-singlet state \af3) = 
(I a) 1 1/3) 2 — \P)i\a)2) /V^- The Hamiltonian is given by. 



n = J2 I dri^ii^) (-^^' - Ma + i?a) V'a(r) + \Y. I dridf2i^i{riWp{T2)Uap^5{ri - r2)V7 (^"2)^5(^1) (26) 



r 



where Ha is the chemical potential of the a-component 
and f/„^^5(f) = /(r)<5a5^/3-y +fl(r)(a|5i|,5) • {(3\S2\l). 
Ea is the energy of hyperfine-Zeeman state \a). Note 
that to the extent that the particle number in any one 
hyperfine-Zeeman level is conserved in the absence of 
the laser beam, that is, if we neglect any decay of 
atoms from one hyperfine-Zeeman state to another, we 
have only two independent chemical potentials, /xi and 
/X2 = /i4, corresponding to the two separately conserved 
quantities A^i = N/2 and iV2 + iV4 = N/2. To address 
the population of the closed channel in a many-body 
system, we look at the equation of motion for a prod- 
uct of two Fermi operators, tpa{ri,t)tpf3{f2,t). 



^^ d{tpa{fi,t)lp0{f2,t)) 

dt 

Vl , El ^2 

2m 2m 



(27) 

M/3 + -£^/3 ) V'a(ri)#(r2) 



+ ^t/a/357(n - f2)V'7(ri)V'5(r2) 



7<5 



+ 



I3'l5-' 

Note that since both /(r) and g{r) are short-range func- 
tions of order ro, it is clear from Eqn.(27) that the 



conversion from an open channel pair state to a closed 
channel molecular state occurs only at short distance, 
i.e., |ri — r2| ^ ^o- It follows then the last two terms 
in Eqn.(27) are of minor importance as compared with 
the other terms since they involve another coordinate 
f which should be close to ri or r2 and thus bring 
up extra factors of kpTQ. They provide either an effec- 
tive background potential or introduce pair states other 
than the ones under consideration (|12) and |14)) which 
are relatively unimportant and thus not of concern here. 
In the following, we shall neglect the last two terms in 
Eqn.(27). Now, taking Eqn.(27) to act on the groimd 
state or thermal ensemble, we find the coupled equation 
of motion of a state with two holes in it. Let us denote 
this state by (f>o{'Pi,'p2) and (/>c (n. , ^^2 ) , where subscript 
o refers to a = 1, /3 = 2 of the open channel and c refers 
to a = 1, = 4 of the closed channel. We assume that 
the rotational degrees of freedom of the closed channel 
molecules are not excited at low temperature and re- 
main in a relative s-wave state. In that case, since W{r) 
is in fact isotropic in space, only the s-wave components 
of (l)o{T\,r2) are important in discussing the population 
of the closed channel molecules. On the other hand, 
only s-wave pair states in the open channel can be con- 
verted by a short-range potential W{r) to closed chan- 
nel molecules, as is clear from the structure of Eq.(28). 
By performing a Fourier transform with respect to the 
center of mass coordinate 2R = r 1 -|- r2 and time, we 
find the following coupled equation. 
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w H h Ml + M2 - -E'^ - K>(?0 (l^o{r; K, w) 

m ' 



w H \- /J-i + IJ-i - Ef> - Sc - eo - Vc{r) ] (j)c{r; K, u) 

m ' 



(28) 



Here — eo is the energy of the molecular state in the 
closed channel relative to its asymptotic energy E2 + E4 

when the two atoms are far away from each other, 
Ej^ = h^K'^/Am is the center of mass kinetic energy of 

a pair of atoms, and Sc = E4 — E2 — eo is the so-called 
bare detuning from the Fcshbach resonance. In case 
of a broad Feshbach resonance, it is much larger than 
the many-body energy scale, in particular, 6c '> cf- 
Even though the form of the coupled equation (28) is 
the same as that for the two-body case, the many-body 
physics does play an important role as it determines the 
normalizations for the function (/>o(n., ^^2) and (t>c{ri, r2)- 
Let us note one feature of Eq.(28) which is conceptually 
important: Since the inter-channel coupling depends 
only on the relative coordinate, so that the center of 
mass momentum K is a good quantum number, we con- 
clude that the pair distributions in the open and closed 
channel are connected by Eqn.(28) due to the super- 
position nature of the open channel pairs and closed 
channel molecules. It is thus in general not permissi- 
ble to assign independent momentum distributions to 
the closed channel molecules and open channel pairs; 
specifying either one of them suffices to fix the other 
through Eqn.(28). Also note that if we neglect, as we 
shall do later, the relatively unimportant factor Ej^ as 
compared with Sc, the coupled equation (28) is identi- 
cal for different i?-states. That implies that, whatever 
the center of mass momentum K is for the open chan- 
nel pair state, the inter-channel coupling always induces 
the same amount of closed channel molecules associated 
with it. The irrelevance of finite .K^-states in discussing 
closed channel molecule formation can again be under- 
stood as a result of its high energy character, namely, 
the process occurs only at short distance, of order ro 
and therefore, many-body physics is quite incapable of 
modifying it. 

Let us then introduce the Green function for the 
closed channel equation. 



where G{r,r') is given by, 



X*n{r)Xn{r') _ x5(r)xo(r') 



ui + eo 



(30) 



where xo('') is the normalized eigenfunction in the 
closed channel with energy — cq. Use Eqn.(30) to in- 



tegrate the closed channel equation, we find, 



,{r;K,ui) = 



(31) 



LO + + fl4, - Ej^ - Sc 

X J dr'xo{r)xo{r')W{r')Mr';K,u). 

Again, we see that since W{r') is short-ranged, the in- 
tegration only picks up the short-range part of the pair 
wave function. An important question is the appropri- 
ate value for w. It is clear that the state obtained by 
removing two particles in states |1) and |2) does not in 
general correspond to the eigenstate(or thermal equilib- 
rium) for the N — 2 particle systcini. It is however, clear 
that u) will be centered around — (/xi + IJ.2), correspond- 
ing to the energy difference between the ground states 
for the N — 2- and A^-particle state. The spread of cu 
will be in general smaller than the Fermi energy even at 
resonance. In the case of a wide resonance, it is known 
that Sc is much larger than the many-body scale, so if 
we approximate the denominator in the above equation 
by Sc and make a Fourier transform with respect to K, 
we find 

Mr,R) = -j J dr'x*o{r)xo{r')W{r')Mr',R)- (32) 

This implies that the number of molecules in the closed 
channel Nq is given by. 



= J drdRMr,RTMr,R) 



(33) 



dRdr'dr"K{r', r")0o(r', i?)</'*(r", R) 



where K{r',r") = Xo{r')Xo{r")W{r')W{r"). It is clear 
that the kernel K{r',r") is a short range function in 
both r' and r", and that 4>o{r' , R)4)*{r" , R) corresponds 
to the s-wave part of the following density matrix, 

(^4{R+-^)^Pl{R-^)^P2{R-<2)MR+^))■ Using 
the same decomposition as before, we find, 

2 



/c 



Nc_ 
N 



kpCKiasM^^T), (34) 



where we have used Eq.(9) and defined Cif(as) = 
/ dr'dr"x{r')*K{r' ,r")x{r"). Noting that the dimen- 
sion ofCK{as) is given by [£']^[Z/], we can define a length 
scale Ic by 



= dr'dr"x{rrK{r',r")xir").{35) 
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Ic is entirely determined by the two-body physics. Wc 
can rewrite the molecular fraction in the closed channel 
as, 



(36) 



At unitarity, fc scales with kp- 

In the following, wc shall illustrate the above general 
considerations in the 'naive' BCS-ansatz, properly gen- 
eralized to include the closed channel component. Note 
that one of the spin state(|l)) is common to the open 
and closed channels, 

\BCS) = Y: {us + vgalal,^ + ^^.a^a^.,,) I^ac)(37) 



where u^, vj: and wj: are the usual variational parame- 
ters, satisfying ju^rp -I- -|- = 1. is the cre- 
ation operator for particle in hyperfine-Zeeman state \i) 

with momentum fc. The corresponding pair wave func- 
tions which are relevant in Eq.(28) are F2 = u^v^ and 

f| = Uj:Wj:. Let us concentrate only on the K = pair 
state in the system since it corresponds to macroscopic 
occupation in the BCS state(cf. also discussion after 
Eq.(28)). Denoting the Fourier transform of F3 and 

FS as F°(r) = EfeWfc^fce'^^-'' and F'^{r) = "feW^fce^'"^ 
and setting u; = — (/ii + ^2), we find that the coupled 
Eqs.(28) take the form. 



Voir) F°{r) = W{f)F''{f) 



m 



m 

- eo - K(r) ) F\f) = W{r)F\r). 



This coupled equation is exactly the same form as that 
in the two-body case(see for example, Ref.[12]). We can 
follow the derivation there or more straightforwardly: 
we can replace ipoir, R) with Q,~^/'^F°{r), where Qr^^"^ 
accounts for the center of mass motion of the pair and 
F°{r) is the radial part of the pair wave function F°{r). 
The spatial dependence of F°{f) is given, within the 
crossover model [38], by 



F^f) 



mA 1 — r/ttg 



(38) 



We find that the density of atoms in closed channel, rig, 
is given by, 



Ir. 



-V 



(39) 



Using n — kp/Sn'^, we obtain the fraction of particles 
in the closed channel fc, 



fc = — = T^kplc — 
n Id Xep 



(40) 



If we compare the above equation with Eq.(36), we 
find /i(C,t) = 37r/16(A/ei;-)2. Thus within the 'naive' 
ansatz, the fraction of particles in the closed chan- 
nel is proportional to and moreover, at resonance. 



scales with kp. In the extreme BEC limit, we know 
A = Aep/ y/Sirkpas and thus. 



f =^ 

Jc — ) 

as 



(41) 



independent of many-body physics, as it should be in- 
tuitively. 

Before ending this section, let us make contact with 
the work in the literature on the problem of the closed 
channel fraction. In the work by Javanainen et.al.[42], 
it is assumed that the Feshbach induced bosons in the 
closed channel are condensed in the K = state. 
This can be regarded as a limiting case of the calcu- 
lation by Chen ei.aL[40] in which Feshbach molecules 
are included in a non-zero temperature generalization 
of the conventional 'naive' ansatz. The conclusions ob- 
tained in [40] are in agreement with our general analysis. 
For example, it is shown in Ref.[40] that the fraction 
of condensed bosons scales with fcp- at unitarity and 
within their approximation, the number of closed chan- 
nel molecules (named Feshbach molecules in Rcf. [40]) 
is proportional to Ag^(our A in Eq.(40) above), while 
the number of non-condensed molecules is proportional 
to A2g(Eq.(9) in [40]). Apg describes non-condensed 
Fermion pairs, which are of course included in the gen- 
eral definition of the function /i(^, r). Thus, it is clear 
that the general structure of the conclusions is the same 
in both approaches. However, as emphasized before, 
as a result of the coupled nature of Eq.(28), the mo- 
mentum distribution of the open channel pair states 
dictates the momentum distribution of the closed chan- 
nel molecules(Feshbach molecules). While this feature 
is explicit at the Hamiltonian level of the two-channel 
model, it is in general not enforced in the actual calcu- 
lations (sec for example Eq.(94)in Ref.[4]). To illustrate 
the point, let us look at the unitarity limit at T = 
where we know that a fraction of the Fermi pairs is not 
condensed; thus the induced closed channel molecules 
associated with them will have nonzero momentum, far 
from being condensed in the K = state. 



IV. CONCLUSIONS 

By exploiting the diluteness of the ultra-cold Fermi 
gas, we have shown that, in considering various physical 
quantities of the system, it is possible to lump all the 
many-body dependence into a single universal function 
/i(^,r). A particular physical quantity may be univer- 
sal, irrespective of microscopic details(e.5., the form of 
the interaction potential), in which case, one should be 
able to express it entirely in terms of /i(^, r), as in the 
case of the average single particle energy of the system, 
Eq.(17). Other physical quantities are not universal and 
there are explicit dependences on the interaction, other 
than that incorporated in the function /i(^,r). How- 
ever, those dependences can be dealt with using only the 
two-body physics. In this case, it is possible to show the 
universal temperature dependence of the physical quan- 
tities. For convenience, let us summarize these two-type 
of behavior in the following. 
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a), Universal dependence on ^ and t. It is under- 
stood that universal here means that all the interaction 
and temperature dependences are captured in one func- 
tion h{^, t). The primary example is the single particle 
energy of the system, Eq.(17). Physical quantities that 
can be directly derived from energy will be in this cate- 
gory as well, for example, the speed of sound c and the 
chemical potential fi in Eq.(23). 

Universal temperature dependence. In this case, 
the physical quantities will have identical tempera- 
ture dependence inherited from /i(^,r). Those physical 
quantities include the interaction energy of the system 
Eq.(12), the average radio- frequency spectroscopic shift 
Eq.(24) and the molecular fraction in the closed channel 
Eq.(36). 

In actual experiment, there is always an external con- 
fining potential which renders the system inhomoge- 
neous. The question of universality is then more del- 
icate. However, the argument given in Sec. II is still 
valid provided the scale over which the confining po- 
tential varies is much larger than the range of the po- 
tential. This is well satisfied in the experiments. The 
universal fimction /i(^(r),r) will depend on position r 
through local Fermi vector kpir)- The temperature de- 
pendence of the physical quantities listed in category /3) 
above will still have universal temperature dependence 
even in a trap. 

This work was supported by the National Science 
Foundation under Grant No. NSF-DMR-03-50842. 

Note added. In the process of writing this paper, we 
become aware of the recent work of Werner eLaL[43], 
where an analysis similar in spirit is carried out for the 
closed channel molecule fraction. Of particular interest 
is their definition of a universal function which is identi- 
cal to our r), see their equation (14). As compared 
with the paper by Werner et.al., we have made an ef- 
fort to connect different physical quantities together and 
emphasized the universal temperature dependences of 
the physical quantities. In addition, we have discussed 
in detail the physical origin of the universal function. 
See however in their paper for a discussion of the tail of 
the momentum distribution(their section 3.2). 



APPENDIX A: EXPANSION OF GROUND 
STATE ENERGY AROUND UNITARITY 

In this appendix, we discuss how to establish the lin- 
ear dependence of the ground state energy e(^) on ^ 
around unitarity. Before we start with the many-body 
problem, it is instructive to look at the two-body prob- 
lem for guidance. For more details, see Rcf.[12] 

Let us consider two atoms of mass m, interacting 
via a central potential V\{r) which can be tuned by 
a parameter A. The wave function for the relative 
motion xxi^") satisfies the following time-independent 
Schodinger equation [12], 

- 2^^XA(r) + VaWxaW = Ex^{r), (Al) 

where mr = m/2 is the reduced mass. Note that the 
normalization of Xxi^) is arbitrary at the moment. Let 



us fix this by requiring that in the region r ro, 



X\{r) = 1 



(A2) 



where ro is the range of the two-body potential V\{r). 
Consider the critical potential V\^, for which ag = oo 

and denote the corresponding radial wave function by 
Xo- Xo = 1 for r Tq. For zero-energy scattering, we 
have, 



d^ 

2mr dr"^ 



Xo(r) + FA,(r)xo(r) =0 



(A3) 



Now, multiplying Eq.(Al) by Xo (setting = in 

the right-hand side as well for zero-energy scattering) 
and multiplying Eq.(A3) by minus xA) we find, using 
Green's theorem and integrating up to tq, 



^^XA(r)-Xo(r)^^ 

/•OO 



or ar 

m r°° 

= dr{Vx{r)-VxAr))xo{r)xxir).{AA) 

Since both Vx{r) and Vx^{r) are short-range functions, 
we can safely replace x\{i~) with xo('') around resonance 
since they are identical for r < rp. Using the explicit 
form of Xa(^) Etnd Xo(^) we find, for infinitesimal change 
of A, 



Sa7^ = - 



dr 



dVxjr) 
dX 



1X0 (r)P 6\A5) 



Now, let us consider the many-body case. We shall 
be interested in a system with N/2 spin up atoms with 
coordinate denoted by Xj and N/2 spin down atoms 
with coordinate denoted by y^, i = 1, 2, • • • , N/2. The 
interaction between spin up atom i and spin down atom 
j takes the form, 



Vx{\xi-yj\ 



(A6) 



where A is a tuning parameter as in the two-body case. 
It determines the asymptotic behavior of the many- 
body wave function in the range ro <C r < a^jkr. We 
denote the corresponding spatial many-body wave func- 
tion as '^\{xi,X2,- ■ ■yi,y2' ■ ■)■ In general, one is not 
allowed to write down a pure spatial wave function with 
spin part totally decoupled from it. However, since we 
are only interested in the energetics of the system, for 
which the spin index is only a spectator, we shall not 
write the spin component explicitly. As in the two-body 
case, the many-body wave function can be normalized 
in such a way that for ro <C \xi — yj\ < ag, kp, 

lim *(xi,X2,---yi,y'2---) (A7) 

ro^\Si-yj\<a,,,kp^ 



1 



1 



47r \Xi - yj\ 



where is a function of variables other than Xi and yj . 
We say that the many-body system is on resonance if 
as OO. as should be regarded as a parameter of the 
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theory which can be tuned in the experiments by the 
external magnetic field. Numerically, the value of as 
defined in Eq.(A7) must be essentially equal to that in 
the two-body case in the same magnetic field, as the 
discussions in Sec. II would imply. 

To put our system in a finite volume such that the 
density of particles n = N/Q is kept constant, we shall 
introduce the characteristic function of the volume O, 

Xn(a?) = l,if X G O, otherwise zero, (A8) 

then we can enforce the condition of constant density 

I 



through an external one-body potential U{x), 

U{x) = Uo{l-Xn{x)) (A9) 

where Uq is a large constant representing the hard wall 

such the many-body wave function vanishes outside the 
region Q and on the boundary of f2, denoted by 



(AlO) 



^'(.T, y) is a short-hand for ^'(iri, X2, • • • j/i, ?72, • ' ' )• The 
time-independent Schodinger equation takes the form. 



V- ^ ft' w ft' w 



In writing the above equation, we have neglected the in- 
teraction potential between parallel spins, correspond- 
ing to the Fock energy in a many-body system. This 
is certainly negligible as compared with the Hartree 
term between anti-parallel spins, which has been in- 
corporated in the above expression. We now follow 
the same recipe developed for the two-body case. We 
write another equation corresponding to A = Aq, i.e. 
corresponding to ttg ^ cxd at resonance. Let us also 
denote the corresponding energy by Eq and the wave 
function by \l/o(a?, y). Then by multiplying each equa- 
tion with ^\ or respectively and subtracting against 
each other, we find. 



E\ — En 



JdxdySVx{xi-yj)^x^o 
J dxdy^x^o 



(A12) 



here 5Vx{xi - yj) = Vx{xi - yj) - Vx^ixi - yj). We 
have used Green's theorem and the fact that the wave 
function vanishes at the boundary of the volume f2. In 
differential form, we have, 



dEx 
dX 



J dxdy 



dVxix,-yj) 
dX 



(A13) 



Now, A determines the scattering length (A) and thus 
if we use ^ = — (fc_FOs)~^, we find, by using Eq.(13), 



dm 

d^ 



da] 



a\ 

J dxdy^x'i'o 

2 



(AM) 



^ _h^j^^E.Jdxdy^\^o\ 



2"^ Io^dr^\Xo(rW JdSdy^x^o 

Notice the resemblance of Eq.(A14) with Eq.(16). In 
fact we have merely managed to express the function 
h{^) in terms of the many-body wave function. We have 



replaced ^x with in the numerator since is a 
short range function. The integral involving in the 
numerator only picks up the short-range contribution 
from the probability distribution function \^Q{x,y)\'^ 
and thus apart from an normalization constant which 
is finite, cancels off the factor dr^^\xo{r)\^ in the 
denominator. Thus, all the sensitive short-range de- 
pendence has disappeared in the above expression and 
we are left with quantities that are independent of 
short-range complications. Note now that the factor 
/ dxdy'^x'^o approaches a constant as ^ ^ since it 
is merely the normalization factor for the wave func- 
tion ^'q. This is why, from a many-body wave function 
point of view, the definition of the function r) is 
universal to the dilute Fermi gas system, irrespective of 
its short-range potential. The complicated expression 
one the R.H.S. of Eq.(A14) reduces to a simple com- 
bination of normalization constants and thus remains 
well-defined for ag ^ 00. We have thus established the 
linear dependence of the energy on ^ around resonance. 

Finally, the extension of the above argument to fi- 
nite temperature is straightforward. In the case of ther- 
mal equilibrium, we instead consider a distribution of 
eigenstates |n) with energy En, each of them occurring 
with probability given by the usual Boltzmann factor 
^-En/ksT ^ One can repeat, word by word, the above 
derivation and thus we can conclude that at finite tem- 
perature, the ^-dependence of energy away from reso- 
nance is linear. This of course, assumes that the tem- 
perature is quite low so that its effects on the short- 
range; wave function is irrelevant. The argument can be 
extended also to a non-equilibrium situation where the 
probability of state \n) is given by p„, X^Pn ^ 1- How- 
ever, it is practically useless since the characterization 
of the system as 'away from the resonance' is ambigu- 
ous and one is not likely to obtain any useful conclusions 
from the argument. 
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